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Abstract 

Using 121.4 fb~^ of data collected with the Belle detector at the T(5S') reso- 
nance at the KEKB asymmetric-energy e^e~ collider, we report evidence for the 
— > A+Att^ decay mode with a measured branching fraction (3.6 ± l.l[stat.] 

+°-^[syst.] ± 0.9[A+] ± 0.7[iV-go]) x 10"" and a significance of 4.4 standard de- 
viations. This is the first evidence for a baryonic B'^ decay. 



1. Introduction 

Results on baryonic B-meson decays obtained by Belle [IHl] and BaBar M- 
[T3 | have increased experimental and theoretical interest in such processes (isj . 
B- meson decay modes with two- [E S 0, HI i three- 0-0, i-im and even four- 
0, 0; [13413 ^^'^ five-body Tl'] final states have been observed. The measured 
branching fractions clearly follow a hierarchy that depends on the final-state 
multiplicity: two-body channels have smaller branching fractions compared to 
multi-body ones. In addition, most three-body baryonic _B-meson decays have 
a near-threshold peak in the invariant baryon-antibaryon mass spectrum. This 
effect was investigated in Ref. [13]. In this paper, we report the first evidence 
for the B^ ^ A+Att" decay and compare the measured branching fraction with 
that for a similar channel, B~ A^pTT~ [4|, where the s-quark of the decay 
under study here is replaced by a u-quark [l5|. 



2. Data Sample and the Belle Detector 

The data for this analysis were taken with the Belle detector at the e~^e~ 
asymmetric-energy collider KEKB [3l at the T(5S') resonance. The integrated 
luminosity of the sample is 121.4 fb~^ and corresponds to (7.1 ± 1.3) x 10^ B'^B^ 
meson pairs [t?] produced in three T{5S) decay channels: B*"B*'^, B*^B'^ and 

The Belle detector is a large-solid-angle magnetic spectrometer that consists 
of a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an ar- 
ray of aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement of 
time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter 
(ECL) comprised of CsI(Tl) crystals located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron fiux return located outside 
the coil is instrumented to detect mesons and identify muons (KLM). The 
detector is described in detail elsewhere [18| . 
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3. Selection Criteria 



We use selection requirements that were previously used for baryonic B- 
meson decay analyses Charged tracks, except those from and A, are 
required to originate within 0.25 cm in the radial direction and within 1 cm 
along the beam direction from the interaction point (IP). We distinguish a 
charged particle of type A from one of type B [A and B being tt, K or p) based 
on likelihood values C{A) and C{B) derived from the TOF, ACC, and dE/dx 
measurements in the CDC. 

mesons (A hyperons) are reconstructed in the — >■ tt+tt^ (A — pn~) 
decay mode by fitting the pion {p and tt) tracks to a common vertex, demand- 
ing an invariant mass in an interval of ±10 MeV/c^ [~ 3cr] (±4 MeV/c^ [w 3ct]) 
around the nominal K'g (A) mass value [19] and applying the following require- 
ments: 

• the distance of closest approach between daughter particles at the decay 
vertex should be less than 3 cm; 

• the distance between the vertex position and IP in the plane transverse 
to the beam direction should be greater than 0.01 cm; 

• the angle a between the Kg (A) momentum vector and the vector pointing 
from the IP to the Kg (A) decay vertex, measured in the plane transverse 
to the beam direction, should satisfy cos a > 0.99; 

• the vertex fit should have < 100 (10) for K^ (A). 

A sample of A+ hyperons is reconstructed in the A+ pK^n^, A+ — ^ pKg, 
and A+ Att^ decay modes. We apply a mass requirement on the recon- 
structed A^ candidates, demanding the invariant mass be within the 10 MeV/c^ 
(« 3(t) interval around the nominal mass value [19l] . 

4. B°-meson Reconstruction 

We fit the A+ and A momentum vectors and the tt track to a common i?^ 
vertex. To reject backgrounds including displaced tracks (e.g., unreconstructed 
Kg and A decays), we impose a loose requirement on the vertex-fit quality for 
the i?^. Signal candidates are identified by two kinematic variables computed in 
the T{5S) rest frame: the beam-energy constrained mass Afbc = \J E'^ca.m ~ P%o 
and the energy difference AE = Ego — i?beamj where i^boam is the beam energy, 
and Ebo and pbo are the energy and momentum, respectively, of the recon- 
structed B^ candidate. For the T(5S') — > B'^B^ production channel, signal 
events correspond to a cluster at (tobo, 0) in the Mhc vs. AE plane. For the 
T(5S') ^ BfW, [BfB*°] channel, the photon from the Bf B^j is not re- 
constructed and so signal events cluster at {{mgo + mB*o) {ttibo — tob*o)/2) 
[(m^.o, m^o — m^.o)] in the Mbc vs. AE plane. We retain B^ meson candi- 
dates with Afbc > 5.3 GeV/c^ and |Ai?| < 0.3 GeV for further analysis. 
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To suppress e+e" cc background, we require that the ratio R2 of the 
second and zeroth Fox- Wolfram moments be less than 0.5. We also specify 
that the angle Othrust between the thrust axis of the candidate in the T{5S) 
frame and the thrust axis of the rest of the event satisfies |cos6thrust| < 0.85. 
The mass window of the A+ candidate, R2 and Othrust requirements are opti- 
mized by maximizing a figure of merit (FOM) Ngig/ -y^^Vbkgd, where iVgig is the 
expected number of signal events from Monte Carlo simulation and A'bkgd is the 
expected number of background events estimated from the AE sidebands in the 
data. 

Signal Monte Carlo samples of 120000 events each for different B^-meson 
production modes and A+ decay channels are used to evaluate the response of 
the detector and determine its efficiency. Events are generated using the EvtGen 



program; the detector response is simulated with GEANT [21[. We model the 
— A+Att" decay according to phase space. 



5. Fit Procedure and Results 

We apply an unbinncd extended maximum likelihood fit simultaneously to 
the three two-dimensional Afbc vs. AE spectra, corresponding to different A^ 
subchannels. Signal and background distributions are parameterized separately 
for all subchannels, taking each function to be the product of shapes in Mbc 
and AE. 

The_contribution of the B° production channel C (C being BfBf, BfB° 
or B'^B^) is parameterized by a two-dimensional Gaussian with parameters 
determined from the Monte Carlo simulation. The typical resolution in Mbc is 
3.6 MeV/c^ and in AE is between 8.8 MeV and 9.9 MeV. The number of signal 
events for the channel C is written as: 

^c'' = /C%0-+A+A7r- ^A+^A7r+ ^A-^prr- '^c' ^ 

where fc is the probability for the B^-meson to be produced through the channel 
C and e is the reconstruction efficiency that is determined from the Monte Carlo 
simulation. For the fractions fc, we use the following values 17[: /g.o-g.o = 



(87.0 ± 1.7)%, /g.oBO = (7-3 ± 1.4)%, and /^o^o = 1 - fs^oBf " /bjobo- 
The — > AjfATT" branching fraction is a common parameter shared among 
subchannels, while the world average values are used for the intermediate 
branching fractions. The average reconstruction efficiency is found to be 12.5% 
for the pKn, 5.9% for the pKg, and 8.7% for the Att subchannel. 

Background shapes are described with an ARGUS threshold function 
in Mbc and a linear function in AE. We exclude the AE < —150 MeV region 
from the fit to avoid contributions from possible B^ A+Att^tt" decays, where 
the 7r° is not reconstructed. This cutoff is verified by Monte Carlo simulation 
of B2 ^ A+A7r-7r° and ^t'^P~ decays. 
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Figure 1: A/bc and AE projections for B° — AjfATr followed by At — >■ pK'n+ (left 
column), — > pK^ (center column) and A^ — > A7r+ (rigiit column). Af^c spectra (top 
row) are for events in the B*°B*'^ signal region (-71 MeV < AE < -23_MeV) and AE 
spectra (bottom row) of the AJAtt" combinations are for events in the B*^B*'-' signal region 
(5.405 GeV/c^ < M^c < 5.427 GeV/c^). The selection requirements and the fit are described 
in the text. 

The fit gives aB°^ A+Att" brandling fraction of (3.6 ± 1.1) x 10"'', which 
corresponds to signal yields of 20.3, 3.0 and 1.9 events for the A+ pif ~7r+, 
A+ — >■ pKg and A+ Att^ subchannels, respectively. The statistical signifi- 
cance of the observed signal is 4.4(T, which is calculated as y^—2\n{Lo/L), where 
Lq and L are the likelihoods with the branching fixed at zero and at the best- 
fit value, respectively. This result provides the first evidence of a baryonic 
decay. Figure [1] shows the one-dimensional Afbc and AE projections for can- 
didates from the BfBf signal region (5.405 GeV/c^ < Mbc < 5.427 GeV/c^, 
-71 MeV < AE < -23 MeV). 

For the systematic error calculation, we change the fixed signal parameters 
of the fit, reconstruction efficiencies, and fractions /a within their uncertainties, 
which gives a contribution of ^i['2%, and change the boundaries of the fit, which 
results in a !'Iq'3% uncertainty. None of these fit variations lower the signal 
significance within the rounding accuracy. We also include a 0.35% per track 
error to account for reconstruction uncertainties, a correlated systematic error of 
2% per p and 1% per tt or to account for the particle identification efficiency, 
an uncertainty of ^g'g % due to the difference between data and Monte Carlo for 
tracks displaced from the IP and errors for all variables that enter into Eqn. ([1]). 
Uncertainties from all sources are summarized in Table [T] and are summed in 
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Table 1: Systematic uncertainties on B{B^l — > AjfAvr ) 



Source 


Relative 
unccrtainty,% 


Fit parameters 


+1.1 

-1.2 


Cutoff 


+0.4 
-0.3 


Tracking efBciency 


±2.1 


Particle identification 


±8.0 


Displaced tracks 


+0.0 
-6.8 




±0.1 




±0.8 


Total 


+8.4 
-12.8 



quadrature. 

Finally, we obtain the branching fraction: 

B{Bl A+Avr") = [i.Q ± l.l[stat.] to:5[syst.] ± 0.9[A+] ± 0.7[Ngo]) x 10"^ 

where the uncertainties due to the A+ absolute branching fractions [l^ and total 
number of iJ^-mesons are shown separately. The B~ A+p7r~ mode, which 
represents a similar decay channel in the Bu-meson sector, has a branching 
fraction of (2.8 ± 0.8) x lO"-* 

To study the observed shapes of the signal and background, we use only the 
A+ pK~Tr^ subchannel as it contains the major portion of the signal. First, 
we examine distributions in the A+ sidebands from 20 MeV/c^ < |M {pK~TT^) — 
m^+l < 50 MeV/c^ and find no peaking structures in the signal area. Monte 
Carlo samples of known T{5S) decays that have six times the statistics of the 
dataset are analysed using the same reconstruction procedure and requirements 
as described above. No hints for the peaking structures in the signal Afbc and 
AE variables are seen. Finally, we check T(5S') -+ B'^*'' B^*\t:) processes [3] 
in which the B^ decays to A+Att". About 120000 signal events are generated 
and then analysed using the same reconstruction procedure and requirements 
as described above. We find no peaking structures in the signal region, while 
a significant background contribution in the /S.E < —200 MeV region is seen. 
From the above checks, we conclude that the signal peak stems indeed from 
B° A+Att^ decays. 

To investigate the possibility of a threshold enhancement, which is com- 
mon in baryonic decays of Bu.d, including the related decay B^ — s> A+p7r~ [3], 
we extract the signal yield in baryon-antibaryon mass bins and, after total re- 
construction efficiency corrections, obtain differential branching fractions as a 
function of M (A+A); these are shown in Fig. [2ja). A fit with a phase-space 
Monte Carlo simulated distribution with floating normalization gives a statisti- 
cal compatibility of 19%. We repeat the same procedure for other two-particle 
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Figure 2: Differential branching fraction as a function of (a) M ^AJT Aj, (b) M ^AJT tt j, (c) 
M (Att-). Points are the data; the solid histogram is the result of the phase-space fit. The 
vertical axis unit is 10~''/(400 MeV/c^). 



invariant masses: A+7r~ [see Fig. [2][b)] and Att^ [see Fig. [H^c)], finding 41% 
and 22% compatibility with the phase-space hypothesis. In general, we observe 
a hint of the typical near-threshold peak in the baryon-antibaryon mass spec- 
trum, while an enhancement to the right could be explained by some decay 
resonant structures. However, the present statistical precision is not sufhcient 
to investigate this effect definitively. 



6. Summary 

In conclusion, we report the first evidence for the — ^ A+Att^ decay 
and measure its branching fraction to be (3.6 ± l.l[stat.] +°;^ [syst.] ± 0.9[A+]± 

0.7[N-^o\j X 10^"* with a 4.4(T significance, including systcmatics. The observed 

— > A+Att^ process represents the first instance of a baryonic decay. 
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